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ABSTRACT 


An  analytic  model  of  long,  propagating  free  wave  perturbations  to 
an  established  upper  ocean  density  front  is  developed.  The  primary 
object  of  the  model  is  to  illuminate  basic  frontal  wave  mechanisms  for 
possible  subsequent  use  in  more  sophisticated,  numerical  models.  The 
model  is  of  the  barotropic  class  but  has  ageostrophic  dynamics  because 
of  the  basic  state  adopted,  essentially  Stommel's  (1966)  model  of  the 
Gulf  Stream  with  uniform  potential  vorticity  and  order  one  Rossby  number. 
The  model  assumes  inviscid  dynamics  apart  from  a  narrow  dissipative  zone 
adjacent  to  the  surface  front.  The  latter  exerts  a  bulk  effect  on  the 
larger  inviscid  zone,  especially  in  generating  small,  but  finite,  cross¬ 
stream  flow  in  the  basic  state.  For  zero  cross-flow  the  resulting  waves 
are  stable,  have  downstream  phase  speeds  that  are  slow  compared  to  the 
current  speed  and  that  increase  with  frequency,  and  have  anomalous 
dispersion.  The  phase  speeds  compare  well  with  the  analysis  of  observati 
of  propagating  Gulf  Stream  meanders  by  Halliwell  and  Mooers  (1983).  For 
finite  cross-flow  the  waves  grow  slowly  in  the  downstream  direction  when 
flow  is  out  of  the  current  and  decay  when  flow  is  into  the  current.  The 
rates  of  growth  or  decay  are  independent  of  wavelength.  Corresponding 
net  growth  or  decay  in  the  wave  kinetic  energy  is  produced  by  action  of 
the  cross-correlation  wave  Reynolds  stress  against  the  lateral  shear  of 
the  basic  state  current. 


1,  Introduction 


Wave  perturbations  to  upper  ocean  density  fronts  are  conmon  mesoscale 
features,  best  known  from  decades  of  observations  of  the  Gulf  Stream 


front  downstream  of  Cape  Hatteras  where  they  have  been  called  meanders. 

Their  major  properties  are  low  subinertial  frequencies  and  long  wavelengths. 
For  the  Gulf  Stream  typical  ratios  of  downstream  wavelength  to  width  of 
the  basic  current  are  about  10  (Halliwell  and  Mooers,  1983).  Perhaps 
most  puzzling  is  the  slowness  of  the  observed  phase  speeds,  typically 
only  a  few  percent  of  the  basic  current  speed.  These  waves  also  show  a 


strong  tendency  toward  downstream  spatial  growth  over  scales  comparable 
to  the  wavelength  itself.  Further  background,  including  an  extended 
discussion  of  previous  models  of  related  wave  dynamics,  is  given  in  my 
recent  paper  on  the  subject,  Garvine  (1983),  hereinafter  termed  683. 

The  model  I  develop  here  is  an  extension  of  G83  and  has  the  following 
seven  major  characteristics: 

1.  The  model  is  effectively  barotropic  in  that  only  the  upper  layer  of 

uniform  buoyancy  is  active;  this  restriction  is  in  contrast  to 
baroclinic  models,  such  as  Orlanski's  (1969)  where  upper  and  lower  ^ 
layers  are  both  active  and  coupled.  ( 

2.  The  undisturbed  or  basic  state  is  a  simple  extension  of  Stommel's 
(1966,  p.  112)  model  of  the  Gulf  Stream.  The  potential  vorticity  P 
is  assumed  to  be  uniform.  Stommel  reasoned  that  since  observations 
of  the  source  water  for  the  Gulf  Stream  in  the  western  Atlantic 
showed  nearly  uniform  P,  its  subsequent  advection  into  the  current 
where  inertial  effects  dominate  should  yield  a  current  of  uniform 
P,  also.  Very  recently  this  reasoning  has  been  significantly 


supported  and  strengthened  by  the  theoretical  work  of  Rhines  and 
Young  (1982  a,b).  Indeed,  maps  of  P  presented  by  McDowell,  et  al . 
(1983)  for  water  above  the  thermocline  and  below  the  direct  reach 
of  atmospheric  influence  show  remarkable  uniformity  over  the  wind 
gyre  of  the  North  Atlantic. 

As  a  consequence  of  (2),  the  model  has  ageostrophic  dynamics, 
because  the  strong  downstream  frontal  jet  of  the  basic  state  has  a 
Rossby  number  of  0(1),  i.e.,  relative  and  planetary  vorticities  are 
of  the  same  order. 

The  upper  layer  flow  is  taken  as  inviscid  apart  from  a  narrow, 
dissipative  zone  adjacent  to  the  surface  front.  There,  however, 
local  turbulence  production  is  likely,  so  that  P  cannot  be  conserved 
along  streamlines.  Johns  and  Watts  (1982)  show  a  vertical  section 
of  P  contours  across  the  Gulf  Stream  downstream  of  Cape  Hatteras 
which  displays  little  change  in  P  above  the  thermocline  except  for 
a  small  area  near  the  surface  front  roughly  coincident  with  cyclonic 
shear  in  the  downstream  current.  Assaf  (1977)  developed  a  steady 
state  model  for  the  Gulf  Stream  where  vertical,  interfacial  friction 
operated  wherever  the  local  densimetric  Froude  number  exceeded  a 
threshold  value.  His  numerical  calculations  showed  realistic 
cyclonic  shear  zone  structure.  Here  the  basic  concept  of  Assaf s 
model  is  employed,  but  not  the  details,  since  they  would  prevent 
analytic  description,  even  in  a  steady  state.  Instead,  as  in  G83, 

I  stipulate  that  wave  perturbations  merely  displace  or  wrinkle  the 
inner  dissipative  zone  without  otherwise  altering  its  structure. 

Its  bulk  effect  on  the  outer,  inviscid  zone  appears  as  constraints 
on  interface  depth  and  cross-stream  velocity  at  the  common  boundary. 


Its  primary  influence  is  the  imposition  of  a  small,  cross-stream 
basic  state  current  upon  the  inviscid  zone  required  to  balance  the 
flux  of  mass  across  the  frontal  interface  in  the  dissipative  zone 
induced  by  turbulent  entrainment. 

5.  The  basic  state  in  general  has  a  small  cross-stream  current  u,  a 
result  of  this  entrainment.  For  u=0  waves  are  stable,  but  for  u>0 
they  grow  downstream,  while  for  u<0  they  decay. 

6.  Only  linear  (small  amplitude)  waves  are  treated. 

7.  The  model  is  developed  for  the  long  wave  limit  only. 

In  G83  I  treated  a  similar  problem.  This  paper,  however,  differs 
in  three  important  ways.  First,  here  I  treat  u^^O,  whereas  in  G83  u=0. 

Second,  here  I  relax  the  stipulation  of  G83  that  shear  in  the  locally 
parallel  velocity  vanish  at  the  boundary  between  the  dissipative  and 
inviscid  zones.  In  G83  this  constraint  permitted  only  stationary  waves, 
that  is,  steady  state  meanders  with  zero  phase  speed.  Here  the  small 
shear  permitted  allows  propagating  waves.  Third,  here  I  treat  only  the 
long  wave  limit  and  thereby  obtain  simple,  analytic  solutions,  whereas 
in  G83  I  treated  arbitrary  wavelength  but  had  to  obtain  solutions  numerically. 

Martin  (1981)  developed  a  numerical  model  with  similar  characteristics 
to  those  above.  In  applications  to  the  Gulf  Stream  front  he  found  slow, 
downstream  wave  propagation  and  sharp  changes  in  rates  of  wave  energy 
growth  with  u  having  fastest  growth  for  u>0.  Pal  dor  (1982)  developed  a 
similar  model  to  G83  having  uniform  P  and  u=0,  but  he  assumed  inviscid 
flow  up  to  and  including  the  surface  front.  Using  a  Rayleigh  integral, 
he  proved  that  this  flow  was  stable  for  small  perturbations  of  any 
wavelength. 


Section  2  presents  the  model  development,  section  3  the  basic 
state,  and  section  4  the  perturbation  equations.  Section  5  then  gives 
the  wave  properties  when  u=0  and  section  6  for  u^^O.  Section  7  concludes 


the  paper. 


2.  Model  development 


Fig.  1  shows  a  schematic  of  the  model  geometry  in  its  basic  (unperturbed) 
state.  Two  incompressible  layers  are  treated,  a  deep  ambient  water 
layer  of  uniform  density  p^,  and  a  light  water  layer  above  the  frontal 
interface  of  uniform  density  P„-Ap.  The  depth  of  the  latter  is  d*  which 
vanishes  at  the  surface  front.  (Asterisks  denote  dimensional  variables.) 

The  source  of  the  light  water  is  called  the  parent  poo.l  where  the  depth 
is  given  as  dg.  Within  the  parent  pool  the  potential  vorticity  is 
uniform  and  given  by  P*=f/dg,  where  f  is  the  Coriolis  parameter. 

The  model  will  be  limited  in  application  to  upper  ocean  density 
fronts  where  the  flow  is  inviscid  apart  from  a  narrow  region  near  the 
surface  front  termed  the  dissipative  zone  (Fig.  1)  where  d*<do*.  The 
dissipative  zone  represents  the  domain  where  locally  high  turbulence 
production  is  expected  from  such  processes  as  mean  vertical  shear  instability 
owing  to  locally  high  densimetric  Froude  numbers,  or,  equivalently,  low 
gradient  Richardson  numbers.  In  this  domain,  in  contrast  to  the  inviscid 
zone,  potential  vorticity  will  not  be  conserved  along  particle  paths. 

Instead,  the  locally  high  vertical  friction  will  generate  cyclonic 
horizontal  shear,  as  in  the  steady  state  models  of  Assaf  (1977),  Kao 
(1980),  and  Garvine  (1980)  for  the  Gulf  Stream  front.  In  addition, 
vertical  volume  flux  by  turbulent  entrainment  across  the  frontal  interface 


Ambient  Water 


Fiaure  1 . 


Schematic  cross-section  of  the  model  front 
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in  the  dissipative  zone  will  force  horizontal  volume  flux  across  the 
inviscid  zone  into  or  out  of  the  parent  pool,  depending  on  the  net 
entrainment  direction  in  the  dissipative  zone. 

Small  amplitude  free  wave  disturbances  will  be  imposed  upon  the 
basic  state.  In  the  inviscid  zone  these  will  be  determined  analytically 
by  solving  the  perturbation  equations  using  an  asymptotic  expansion  for 
small  wavenumber.  In  the  dissipative  zone  no  explicit  solution  will  be 
attempted  for  the  perturbations,  because  even  the  basic  state  there  is 
too  complicated  for  analytic  solution.  Instead,  as  in  G83,  I  will 
stipulate  that  wave  perturbations  will  merely  displace  the  dissipative 
zone  without  otherwise  altering  its  structure.  Its  bulk  effect  on  the 
larger  inviscid  zone  will  nevertheless,  be  accounted  for  by  inner  boundary 
conditions  on  the  interface  depth  and  cross-stream  velocity  at  the 
common  boundary,  termed  the  inner  boundary  (Fig.  1).  The  primary  effect 
of  the  dissipative  zone  wil"  be  the  generation  of  a  source  or  sink  of 
cross-stream  volume  flux,  which,  in  turn,  will  affect  fundamentally  the 
stability  of  the  waves  in  the  whole  domain. 

Standard  oceanographic  model  restrictions  are  imposed.  The  Boussinesq 
approximation  is  assumed  valid,  wind  stress  is  neglected,  vertical 
gradients  in  horizontal  current  in  the  inviscid  zone  are  absent,  and 
beta  effects  are  neglected.  The  upper  layer  is  assumed  shallow  compared 
to  the  water  depth,  so  that  the  upper  layer  pressure  field  is  in  isostatic 
equilibrium.  The  ambient  water  is  taken  to  be  at  rest  and  the  front  to 
have  zero  mean  velocity  of  propagation  normal  to  itself;  these  latter 
restrictions  may  be  removed  following  G83,  but  have  no  fundamental 
effect  on  the  properties  of  the  frontal  waves  or  their  stability. 
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The  independent  variables  are  the  time  t*,  the  cross-stream  distance 
X*  (Fig.  1),  and  the  downstream  distance  y*;  u*  and  v*  are  the  corresponding 
velocity  components  and  the  potential  vorticity  is  P*.  Following  G83  I 
introduce  scaled  variables  as 


(x,y)  e  (x*,y*)/\,  t=ft* 

(u,v)  =  (u*,v*)/c^,  d=d*/djj* 

P=P*djj*/f 

where  c.=(g'd.*)^  with  g'=g  Ap/p  ,  the  reduced  gravity,  so  that  c.  is 

ID®®  I 

the  linear  internal  wave  velocity  and  xec^/f,  the  internal  Rossby  radius. 
As  these  scaled  variable  imply,  c^  is  the  scale  for  the  current,  x  for 
the  cross-stream  distance,  and  f/d^^*  for  the  potential  vorticity. 

In  these  scaled  variables  the  governing  equations  of  potential 
vorticity,  cross-stream,  and  downstream  momentum  are  respectively  (G83): 


0 

Ot 

V 

Dt  3x 

!5t  ^  ^  ^  W 


(1) 

=  0, 

(2) 

=  0, 

(3) 

where  0/Dt  =  3/3t  +  u3/3x  +  v  3/3y,  the  material  derivative,  and  P  = 
(3v/3x  -  3u/3y  +l)/d.  Here  (1)  states  the  conservation  of  potential 
vorticity  along  particle  trajectories  and  will  be  used  as  an  alternate 
to  the  continuity  equation. 


3.  The  basic  state 


The  flow  state  is  decomposed  in  standard  fashion  into  a  basic  state 
that  represents  the  steady,  undisturbed  flow  and  a  perturbation,  here 
assumed  small.  Thus, 

u  =  u  +  u' ,  V  =  V  +  v* ,  (4) 

d=a+d',  P=P+P', 

where  overbars  denote  basic  state  and  primes  perturbation  variables. 

As  in  G83,  the  essential  model  physics  arises  from  the  assumption 
of  uniform  potential  vorticity  for  the  basic  state  in  the  parent  pool 
and,  thus,  for  the  entire  inviscid  zone  as  well,  an  idea  Stommel  (1966) 
was  first  to  postulate.  Consequently,  F  =  1.  Since  I  further  assume 
the  front  is  straight  and  parallel  to  y  in  the  basic  state,  9/3y  «  3/3x. 
As  shown  in  G83,  the  solutions  for  3  and  v  for  the  inviscid  zone  without 
cross-stream  flow  (u=0)  are  then 

a  =  l-Voe“^,  V  =  1-a  =  Voe'^,  (5) 

where  Vo  =  1-do  with  do  the  scaled  depth  at  the  inner  boundary.  Equations 
(5)  are  the  same  as  Stormiel  (1966,  p.  Ill)  obtained,  except  here  Vo£l , 
i.e.,  the  frontal  interface  need  not  surface  at  the  inner  boundary 

(do>0). 

Small,  but  finite,  basic  state  cross-stream  volume  flux  in  the 
inviscid  zone  will  be  generated  by  mass  entrainment  in  the  dissipative 
zone.  Within  the  inviscid  zone  this  volume  flux  is  uniform,  so  that  in 
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the  absence  of  y-derivatives  we  have  from  mass  continuity  u*  3*  = 
where  is  the  (dimensional)  cross-stream  flow  into  the  parent  pool. 


For  convenience,  let  this  be  represented  by 


^i  'e  ' 


where  denotes  the  sign  of  the  entrainment  (S^  =  +1  for  entra  “nt 


into  the  upper  layer  and  =  -1  for  entrainment  into  the  ambie  layer) 


while  F  =  |u^*|/c^.  is  the  Froude  number  indicating  the  relative  ,  _-ed  of 
the  outflow  (or  inflow).  In  contrast  to  small-scale  density  fronts 


where  F  =  0(1)  (Garvine,  1974;  Kao,  et  al.,  1977),  for  the  large  scale 


fronts  of  interest  here  F«1 ,  typically  10'  .  With  this  nomenclature  we 


then  have 


u  =  Sg  F/a. 


It  may  be  shown  that  u  alters  the  basic  state  solution  of  (1)  -  (3)  only 


at  O(F^).  Because  I  will  include  only  terms  through  0(F)  subsequently, 
(5)  will  still  be  valid  for  3  and  v. 


F  will  emerge  as  the  small  parameter  associated  with  frontal  wave 


stability.  The  basic  state  of  G83,  in  contrast,  assumed  u=0  (or  F=0) 


and  was  associated  with  only  stable  waves. 


4.  The  perturbation  equations 


Substitution  of  the  decomposition  forms  (4)  into  the  governing 


equations  (1)  -  (3)  and  retention  of  only  first  order  perturbation  terms 


^ 
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yields  the  inviscid  zone  linear  perturbation  equations.  The  perturbation 
potential  vorticity  equation  is  especially  simple,  since  (1)  dictates 
conservation  of  P  following  the  fluid,  while  P  =  1  by  postulate.  Thus, 
P'  =  0  everywhere  in  the  inviscid  zone,  yielding: 


1^' 

3x 


ly.' 

3y 


-  d'  =  0. 


(7) 


The  perturbation  x-  and  y-  momentum  equations  are 


3t 


-  v‘ 


S  F  — 
e  3x 


=  0, 


at 


+  V 

3y 


M' 

3y 


(8) 

(9) 


Note  the  appearance  in  (8)  and  (9)  of  both  downstream  advection  terms, 
with  v(x)  as  a  factor,  as  well  as  cross-stream  advection  terms  containing 
SgF/3  =  u;  these  latter  terms  bring  the  small  parameter  F  into  the 
problem. 

Following  G83,  I  adopt  the  wave  forms: 

u'  (x,y,t)  =  -i  [Xo(x)  +  iFX^(x)]e^^,  (10a) 

v'  (x,y,t)  =  [Yo(x)  +  iFY^(x)]e^^,  (10b) 

d'  (x,y,t)  =  [Zo(x)  +  iFZi(x)]e’®,  (10c) 


where  i  =  /-I  and  0  is  the  phase  function  given  by 


0  (y.t)  =  (k-iFK)y--..;t. 


(lOd) 


In  evaluating  any  of  the  perturbation  variables  the  real  part  must  be 
taken.  Here  k  is  the  scaled  wavenumber,  taken  as  real,  while  u  is  the 
scaled  wave  frequency,  shown  later  to  be  real,  and  <  anticipates  the 
possibility  of  wave  instability  (<>0)  or  decay  (<<0)  when  the  wavenumber 
has  an  imaginary  part  for  small,  but  finite,  F.  The  wave  forms  anticipate 
at  u'  will  be  in  quadrature  phase  with  v'  and  d'  when  F=0,  while  a 
primary  effect  of  F>0  will  be  to  change  the  relative  phases  of  the 
variables.  The  eigenfunctions  Xo(x),  X^(x),  Yo(x),  etc,  will  all  turn 
out  to  be  real.  Subscripts  "o"  denote  the  eigenfunctions  for  the 
limiting  case  of  F=0  (no  cross-flow  in  the  basic  state)  while  subscripts 
"1"  denote  corrections  for  F  small  but  finite. 

Substitution  of  the  wave  forms  (10)  into  the  perturbation  equations 
(7)  -  (9)  and  retention  of  only  0(1)  terms  yields  the  following  zeroth 
order  system  of  ordinary  differential  equations. 


o  ■  k  X  -  7  =  n 

K  Ao  Lo  u, 

(k  V  -  C.)  Xo  -  Yo  +  ^  =  0, 

(k  V  -  0,)  Yo  -  (1+^)  Xo  +  k  Zo  = 


(lla) 

(lib) 

(lie) 


These  are  identical  in  form  to  (14)  of  G83  with  w=0,  as  expected,  since 
there  u=0. 

Retaining  terms  through  0(F)  yields  the  first  order  system. 


-  kX^  -  =  -kXo 


(kv  -  a))X^  -  +  37^  =  <vXo  +  Sg 

(kv  -  a))Yi  -  kZ-j  =  k(vYo  +  Zo)  +  7^  . 

d 


(12a) 


(12b) 


(12c) 


•*—  ■* _ 


The  left  sides  of  these  are  luentical  to  (11),  but  inhomogenous  or 
forcing  terms  now  appear  on  the  right  generated  by  the  coupling  of 
finite  cross-flow  and  downstream  wave  amplitude  growth  or  decay  to  the 
F=0  wave  field. 

5.  Long  waves  without  cross-stream  flow 

I  treat  in  this  section  long  waves  without  cross-stream  flow  in  the 
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basic  state,  i.e.,  in  the  limits  k  «1  and  F=0  (or  u=0).  The  system  to 
be  solved  is  (11)  for  k  «1  subject  to  the  boundary  conditions  for 
frontal  "trappea"  waves 

Xo,  Yo,  Zo  ^  0.  for  x-x»,  (13) 

and  to  those  at  the  displaced  inner  boundary  where  x  =  h(y,t).  The 
inner  boundary  displacement  may  be  written,  consistent  with  the  wave 
forms  (10),  as 

h(y,  t)  =  He^’®, 

where  H  is  the  frontal  displacement  amplitude.  The  inner  boundary 
conditions  I  impose,  as  in  G83,  are  that  the  interface  depth  at  the 
displaced  inner  boundary  is  unchanged  from  that  of  the  basic  state,  do, 
while  the  fluid  velocity  component  locally  normal  to  the  boundary  is 
likewise  unchanged.  For  u=0  these  are  (G83),  respectively: 

Zo(0)  =  -VoH  =  -1.  (14a) 

and  Xo(0)  =  -kVoH  +  ojH  =  -k  +  u/Vo,  (14b) 

where  I  equate  VoH  to  unity  for  convenience,  since  the  particular  amplitude 
of  a  free  wave  may  be  chosen  arbitrarily. 


The  mathematical  problem  then  is  to  solve  the  system  (11)  subject 

to  the  boundary  conditions  (13)  and  (14).  This  is  nearly  the  same 

problem  as  in  G83,  but  with  two  important  differences.  First,  the 

stipulation  of  G83  that  there  be  no  shear  of  the  locally  parallel  velocity 

component  at  the  inner  boundary  is  relaxed,  i.e.,  now  only  the  normal 

component  must  match  that  of  the  (unchanged)  dissipative  zone.  Second, 

an  asymptotic  expansion  is  employed  yielding  analytic  solutions  through 
2  2 

0(k  )  for  k  «1 .  The  first  difference  permits  the  existence  of  propagating, 
as  opposed  to  stationary,  waves,  while  the  second  provides  simple, 
analytic  solutions  in  a  limit  of  primary  physical  interest. 

If  one  stipulates  that  the  dissipative  zone  is  merely  displaced  by 
the  wave  perturbations,  but  is  otherwise  unaltered,  as  here  and  in  G83, 
and  also  stipulates  the  no  shear  condition,  as  in  G83,  then  one  has 
Yo{0)  =  VoH  =  1.  (15) 

If  this  is  inserted,  together  with  the  inner  boundary  conditions  (14), 
into  the  y-momentum  equation  (11c)  and  evaluated  at  x=0,  one  finds 
a»=0  for  all  k,  i.e.,  only  stationary  or  steady  state  waves  are  possible, 
as  in  G83.  Conversely,  relaxing  the  no-shear  stipulation,  as  here, 
permits  oifO  or  propagating  waves. 

Appropriate  asymptotic  expansions  for  the  long  wave  limit  are 

Xo(x)  =  kCRo^°^(x)  +  k^Ro^^hx)  +  •••] 

Yo(x)  =  Yo^^^x)  +  k^Yo^^^x)  +  •••  •  (16) 

Zo(x)  =  Zo^°\x)  +  k2Zo^^^(x)  +  ... 

2 

Here  k  is  the  small  expansion  parameter,  while  all  eigenfunctions  on 

the  right,  including  Ro^^^  and  Ro^^^,  are  0(1).  These  forms  thus  anticipate 
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that  Xo  (or  u')  will  be  0(k),  while  Yo  (or  v')  and  Zo  (or  d')  will  be 

0(1). 

Substituting  (16)  into  (11),  assuming  that  u)  =  0(k"')  with  m>1 ,  and 

2 

retaining  only  the  zeroth  order  terms  in  k  yields 


o 

II 

o 

o 

1 

X 

o 

0 

>“ 

“O 

(17a) 

-Yo^°^  +  dZo^°Vdx  =  0, 

(17b) 

vYo^°^-(l+^)Ro^°^  +  Zo^°^  =  0. 

{17c) 

Equation  (17a)  indicates  that  potential  vorticity  changes  are 
absent  as  a  result  of  a  balance  between  the  perturbations  in  lateral 
shear  and  interface  depth,  while  (17b)  indicates  the  cross-stream  momentum 
balance  to  this  order  is  geostrophic.  This  pair  may  be  solved  easily 
and  the  boundary  conditions  (13)  and  (14a)  imposed  to  give 

Yo^°^(x)  =  -Z^^hx)  =  e‘^.  (18) 


Thus,  at  this  order  of  approximation,  Yo(0)  satisfies  (15),  indicating 
zero  shear  of  the  parallel  velocity  across  the  inner  boundary. 

The  downstream  momentum  equation  (17c)  may  now  be  solved  easily  to 

yield 


(19) 


This  identically  satisfies  the  boundary  condition  (14b)  to  zeroth  order. 
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The  simplicity  of  (17c)  and  its  solution  (19),  however,  masks  a 

complicated  downstream  momentum  balance  between  four  of  the  six  terms  of 

the  perturbation  equation  (9):  downstream  advection  by  the  basic  state 

(v  3v'/3y).  the  Coriolis  term  (u'),  cross-stream  advection  by  the  perturbations 

(u'dv/dx),  and  the  perturbation  pressure  gradient  (3d'/3y).  Consequently, 

(17c)  is  highly  ageostrophic,  reflecting  the  0(1)  Rossby  number  of  the 

basic  state,  and  lacks  only  the  local  acceleration  term  and  cross-stream 

advection  by  the  basic  state  of  (9). 

2 

The  first  order  equations  in  k  are 


dYo^^Vdx  - 

+  dZo^^  Vdx  =  -vRo^°^ 
vYo^^^  -(iARo^^^  +  Zo^^^  =  ((u/k^)Yo(°^ 


(20a) 

(20b) 

(20c) 


This  system  differs  from  the  zeroth  order  one  in  three  ways:  the  cross¬ 
stream  flow  now  participates  in  the  vorticity  balance  of  (20a),  downstream 
advection  of  momentum  by  the  basic  state  renders  even  the  cross-stream 
momentum  balance  ageostrophic  in  (20b),  and  the  local  acceleration  is 
added  to  the  downstream  momentum  balance  in  (20c).  The  latter  represents 
the  only  time  dependent  term  in  the  governing  equations  to  first  order 
and  brings  u  into  the  problem.  Consistency  requires  u  =  0(k  )  or  m=3. 

Equations  (20a)  and  (20b)  may  be  solved  for  Yo^^^  and  Zo^^^  subject 
to  (13)  and  the  boundary  condition  Zo^^^(O)  =  0,  consistent  with  (14a), 
to  yield 


,<"(.)  =  4  ixe-\ 


(21a) 
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Vo^^^x)  =  (^  -  I  Vo)e‘^  +  y  Voe"^^  -  ^  xe"^.  (21b) 

If  (20c)  is  now  solved  for  and  evaluated  at  x=0  one  finds, 

Ro^^^(O)  =  U  (3-2vo)  -  ^  . 

dok'^ 

2 

However,  the  boundary  condition  (14b),  if  evaluated  to  0(k  ),  gives 


Ro^^^O)  =  . 

Vok-^ 

Equating  these  gives  the  dispersion  relation 

2  , 

0)  =  ^  (3-2Vo)k'^. 


(22) 


Since  k  is  small  and  cu  =  0(k  ),  these  waves  have  low  subinertial 
frequencies  and  slow  phase  speeds.  Because  the  waves  are  so  slow,  their 
corresponding  downstream  fluid  acceleration  av'/st  enters  the  downstream 
momentum  balance  only  at  0(k  ),  so  that  the  solutions  to  this  order,  as 
well  as  to  0(1),  must  be  found  in  order  to  determine  w. 

For  the  limit  of  vanishing  dissipative  zone  (do-*0  or  Vo-»-l) 
uj->k  /6.  This  limit  agrees  with  Paldor's  (1982)  analysis  of  his  model 
for  long  waves,  as  expected,  since  Paldor's  model  assumed  inviscid  flow 
everywhere  (do=0).  For  small  do  the  more  general  result  of  (22)  is 
quite  close  to  the  limiting  value;  for  example,  when  do=.l  (or  Vo=.9), 
(22)  gives  a)=.162k  ,  a  difference  of  less  than  3%  from  the  limiting 
case. 
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The  phase  and  group  velocities  derived  from  (22)  are 
V  2  2 

c  *  u/k  =  {3-2vo)  k^,  (23a) 

and  Cg  =  3(jL)/3k  =  3c.  (23b) 

Since  Cg>c,  these  waves  have  anamalous  dispersion.  They  propagate  only 
downstream  with  phase  speeds  of  0(k  ),  i.e.,  much  slower  than  the  downstream 
current  speed.  Martin  (1981)  developed  a  numerical  model  of  frontal 
waves  which  featured  a  simple,  active,  upper  layer  and  a  dissipative 
zone  that  also  was  advected  by  frontal  displacements,  but  otherwise 
unchanged,  as  here.  His  basic  state  structure  was  also  similar.  He 
imposed  a  half  wave  sinusoidal  displacement  to  the  basic  state,  simulating 

a  trough,  and  computed  the  subsequent  evolution  of  the  flow.  The  particular 

2 

case  most  related  to  the  present  model  had  an  equivalent  value  of  k  -0.147, 
within  the  long  wave  limit,  kH=0.105  (small  amplitude),  F*0,  and  Vo=2/3. 

He  found  downstream  propagating  waves  of  normalized  phase  speed  c®0.0334, 
while  (23a)  gives  c*0.0181,  somewhat  slower.  He  also  found  that  shorter 
waves  developed  downstream  of  the  initial  wave  which  had  higher  phase 
speeds,  implying  anomalous  dispersion  behavior,  as  in  the  present  model. 

The  low  speed  of  long  waves  (meanders)  of  the  Gulf  Stream  compared 
to  the  swiftness  of  the  current  has  been  often  remarked  (Hansen,  1970; 
Halliwell  and  Mooers,  1983).  Comparison  of  present  results  with  observations 
of  deep  water  Gulf  Stream  meander  properties  is  therefore  of  interest, 
though  such  comparisons  must  be  made  cautiously,  because  the  ratio  of 
upper  layer  to  total  depth  is  not  very  small,  about  1/5. 


The  most  comprehensive  observational  data  set  for  deep  water  Gulf 
Stream  meanders  is  that  of  Halliwell  and  Mooers  (1983).  They  analyzed 
meander  properties  from  a  four-year  collection  of  weekly  charts  derived 
from  satellite  infra-red  images  for  the  region  up  to  1000  km  downstream 
of  Cape  Hatteras  and  found  two  dominant  modes:  very  low  frequency 
standing  waves  with  periods  of  months  to  years  and  downstream  propagating 
waves  with  periods  of  several  weeks.  The  present  results  pertain  to  the 
latter,  propagating  mode. 

First,  I  compare  results  for  one  particular  case,  then  I  compare 
phase  speeds  and  wavelengths  over  a  wide  frequency  range.  To  compute 
dimensional  values  from  the  present  results  for  comparison  the  variables 

_C  _1 

must  be  unsealed.  I  selected  as  typical  of  the  Gulf  Stream  f=9xl0  s  (39 

and  c^=2.5  m/s,  corresponding  to  a  mean  density  difference  of  Ap/p^ 

-4 

=  7x10  and  a  pycnocline  depth  d|j*®900  m;  consequently  the  Rossby 
radius  is  A*27.8  km.  For  simplicity  I  took  do=0  (vo-l).  The  model  then 
gives 

c*  =  1.87  c.  (.Vf)^^^,  (24a) 

A  =  1.87  A(n/f)"^/^,  (24b) 

where  c*  is  the  dimensional  phase  velocity,  A  the  dimensional  wavelength, 

and  n  h  1/T*,  the  reciprocal  period.  For  fi=8  cycles/year  (cpy),  Halliwell 

and  Mooers  (1983)  found  for  the  most  energetic  waves  values  of  c*=9.1 

cm/s  and  .A=360  km;  the  model  gives  c*=9.3  cm/s  and  A=369  km,  corresponding 
2 

to  c=.0373  and  k  =0.224.  They  roughly  estimated  the  group  velocity  c^* 
for  the  frequency  band  4  to  10  cpy  at  17  cm/s,  while  the  model  gives  28 
cm/s.  Results  for  c*  and  A  are  compared  in  Fig.  2.  The  model  gives 


Figure  2.  Comparison  of  the  observational  results  from  Halliwell 
and  Mooers  (1983)  and  the  present  model  for  wave  phase 
speed  c*  and  wavelength  A  for  wave  frequency  Q  up  to 
10  cycles/year.  Circles  denote  the  values  for  the  most 
energetic  waves  and  the  brackets  the  associated  95% 
confidence  intervals  from  Halliwell  and  Mooers. 


2/3 

c*a  Q  and  agrees  well  with  both  the  trend  and  the  particular  values 

of  the  observations.  Alternatively,  one  may  compare  values  of  A, 

although  these  are  not  independent  of  c*,  since  Values  for  olO 

2 

cpy  were  not  compared  because  k  >0.26  then,  so  that  the  long  wave  limit 
is  dubious. 

In  essence,  the  observations  and  model  results  have  the  following 
common  features:  comparable,  slow  phase  speeds  for  given  frequency, 
increasing  phase  speeds  with  increasing  frequency,  and  anomalous  dispersion 
(c*  >  c*).  This  agreement  suggests  that  the  model  provides  a  simple, 
"lowest  order"  description  of  such  waves  in  the  same  sense  that  Stommel's 
(1966)  model,  used  here  for  the  basic  state,  does  for  the  mean  Gulf 
Stream  downstream  of  Cape  Hatteras. 

2 

The  scaled  eigenfunctions  Xo,  Y©,  and  Zo,  computed  through  0(k  ), 

as  well  as  their  0(1)  parts  are  plotted  x  in  Fig.  3  for  the  particular 

2  2 

case  above  where  do=0  and  k  =0.224.  Note  that  the  0(k  )  corrections  are 

2 

modest  even  for  this  relatively  large  k  .  The  trapped  character  of  the 
waves  is  evident  in  the  clear  exponential  decay  for  large  x.  One  may 
more  readily  interpret  the  eigenfunctions  by  taking  the  real  parts  of 
the  wave  forms  (10)  to  give,  for  F=0, 

u'  =  Xo  sin  0,  v'  =  Yo  cos  0,  d'  =  Zo  cos  0,  h  =  H  cos  0, 
where  0(y,t)  is  the  phase  function  given  by  (lOd).  Since  Yo^O,  the 
downstream  wave  velocity  v'  is  in  phase  with  the  frontal  displacement 
h,  so  that  at  wave  troughs  (h=H  or  maximum  frontal  displacement  toward 
the  parent  pool)  the  downstream  flow  is  faster  than  the  basic  state, 
while  at  wave  crests  (h=-H)  it  is  slower.  Conversely,  since  Zo^O,  the 
interface  displacement  d'  is  out  of  phase  with  both  h  and  v',  so  that 
the  interface  shoals  at  the  wave  troughs  but  deepens  at  the  crests.  In 


contrast  u  is  in  quadrature  phase  with  h,  v  ,  and  d  ,  having  maximum 
values  where  h=0  and  zero  values  at  troughs  and  crests.  The  quadrature 
phase  with  v'  means  that  the  Reynolds  stress,  the  mean  of  -u'v'  over  a 
period,  will  vanish,  consistent  with  stable  waves. 


The  velocity  shear  across  the  inner  boundary  may  be  expressed, 
following  (14)  and  (20c),  as 


k^Yo^^^(O)  cos  0  =  k^cos  0  oo/k^Vo  =  c  cos  0/Vo. 


If  this  shear  were  required  to  vanish  for  all  y  and  t,  as  in  G83,  then 

(jj=0=c,  and  the  waves  would  be  steady  and  stationary,  as  in  G83.  There 

2 

the  long  wave  limit,  k  ->0,  corresponded  to  only  a  trivial  solution  for 
the  zeroth  wave  mode.  All  nontrivial,  higher  wave  modes  had  kSo(l), 
so  that  no  long  waves  were  possible.  The  present  results  may  therefore 
be  viewed  as  an  extension  of  G83  for  the  zeroth  mode  or  long  wave  limit 
obtained  by  relaxing  the  condition  of  no  shear  across  the  inner  boundary. 


6.  Lonq  waves  with  slow  cross-stream  current 


I  treat  in  this  section  long  waves  with  slow  cross-stream  current 
in  the  basic  state,  i.e.,  k  «1  with  0<F<<1  (or  u^O).  The  wave  forms 
are  given  by  (10)  where  the  eigenfunctions  to  be  found  now  are  Xi(x), 
Yi(x),  and  Z^(x),  while  the  governing  equations  are  given  by  (12).  The 
boundary  conditions  are  that  the  eigenfunctions  vanish  for  x-«°,  as  in 
(13)  for  the  F=0  flow,  and  that  the  interface  depth  and  locally  normal 
fluid  velocity  at  the  displaced  inner  boundary  match  that  of  ^‘■e  .  '‘jace 
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dissipative  zone.  The  latter  conditions  are  extensions  of  (14)  and  may 
be  written 

Z^(0)  =  0.  (25a) 

Xi(0)  =  K  +  Sg/do^.  (25b) 

Equation  (25a)  follows  from  full  satisfaction  of  the  interface  boundary 
condition  (14a)  by  Zo(0)  itself,  while  (25b)  follows  from  the  additional 
physics  brought  by  downstream  amplitude  change  and  basic  state  cross- 
fl  ow. 

The  boundary  condition  (25b)  implies  X^(x)  =  0(1),  while  (12a, b) 
imply  that  Y-j  and  Z-j  are  0(k).  Thus,  expanding  X-j  as 

Xi(x)  =  x/°^x)  +  x/^^(x)  +  ..., 

and  inserting  it  in  (12c),  the  downstream  momentum  equation,  gives  to 
0(1)  in  k^: 


-(1  +  ^)x/°^  =  <(vYo^°^  +  Zo^°^)  +  ^ 
After  substitution  for  the  F=0  variables,  one  finds 

Q 


dYo 

dx 


(0) 


(26) 


This  satisfies  the  boundary  condition  (25b)  identically. 


Solutions  for  Y-i  and  may  now  be  obtained  to  lowest  order  in  k 
from  (12  a,b).  With 

=  k  S(x),  =  k  T(x), 

one  may  substitute  for  X-j  from  (26)  and  eliminate  S  in  favor  of  T  to 
obtain 


2 

^  -  T  =  2Ke‘^  (2v  +  1).  (27) 

dx^ 

Using  the  boundary  condition  T(0)  =  0  from  (25a),  one  may  solve  this  to 
give 


T(x)  =  <e"^  V  -  ^  Vo  -  x). 


(28) 


which,  in  turn,  yields, 


S(x)  =  Ke"^  (|  Vo  -  I  V  -  1  +  x)  -  —  e"^. 

j  j  3 


(29) 


To  close  the  0(F)  system  for  <,  one  must  continue  the  solutions  to 


O(k^).  Evaluating  the  downstream  momentum  equation  (12c)  to  this  order 
gives 


a  =  V  S  +  T  -  <vY 


(1) 


(1)  ^e  dY^ll 


dx 


Here  the  right  side  is  now  completely  known,  except  for  <.  Evaluating 


this  at  x=0  and  recognizing  that  '(0)  =  0,  since  X-j  ^  ^(0)  already 
satisfies  (25b),  gives 


K  = 


Vo(3-2vo)' 


(30) 


This  is  the  principal  result  of  this  section  and  shows  three  primary 
features  of  the  model  wave  growth.  First,  the  sign  of  <  is  the  same  as 
that  of  Sg.  Thus,  for  5^=1,  when  there  is  basic  state  cross-flow  away 
from  the  front  (upward  entrainment),  the  waves  grow  in  the  downstream 
direction,  indicating  instability;  conversely  for  Sg=-1 ,  when  the  cross- 
flow  is  toward  the  front  (downward  entrainment),  the  waves  decay,  indicating 
stability.  Second,  <  is  independent  of  wavenumber,  so  long  as  the  waves 
are  long;  thus,  there  is  no  preferential  or  fastest  growing  unstable 
wave.  Third,  k=0(1),  so  that  its  dimensional  counterpart  form  (lOd)  is 
<*  =  Fk/X, 

much  smaller  than  i.e.,  the  downstream  growth  scale  will  be  0(\/F), 
many  Rossby  radii. 

The  best  documented  observations  of  growing  frontal  waves  are  those 

of  the  Gulf  Stream  downstream  of  Cape  Hatteras.  Hansen  (1970)  found  <* 

=  (2±l)xl0'^  km“^  and  Watts  and  Johns  (1982)  about  2.5x10”^  km”\ 

Halliwell  and  Wooers  (1983)  analyzed  the  greatest  amount  of  data  for 

propagating,  growing  waves.  For  the  frequency  band  4<n<8.5  cpy  they 

-3  -1 

found  K*  =  (3.2±1.3)xl0  km  with  no  perceptible  variation  with  n. 

The  model  results  similarly  show  no  change  with  a.  For  Vo=l  and  5^=1, 

(30)  gives  k=2.  Thus,  for  x=27.8  km,  as  in  section  5,  <*  would  match 
that  of  Halliwell  and  Wooers  for  F=0.044,  corresponding  to  a  cross¬ 
stream  flow  into  the  parent  pool  of  Uu*  =  Fc-  =11  cm/s,  the  order  of 
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magnitude  one  would  expect.  This  value  of  <*  corresponds  to  a  downstream 
e-folding  distance  of  310  km,  roughly  comparable  to  the  wavelength  A 
predicted  (Fig.  2)  for  this  frequency  band. 

Further  comprehension  of  the  nature  of  wave  amplitude  growth  or 
decay  in  the  model  is  available  through  analysis  of  the  wave  kinetic 
energy.  One  may  derive  an  equation  describing  this  by  first  adding  the 
product  of  (8)  with  3u'  to  the  product  of  (9)  with  dv‘,  then  averaging 
the  results  over  a  wave  period,  and  finally  integrating  over  the  entire 
in viscid  domain.  One  finds 


dKVdt  =  .  A,  .  Wp. 


(31) 


where  K' 


the  total  wave  or  fluctuation  kinetic  energy. 


the  rate  of  working  by  wave  crosscorrelation  Reynolds  stress  against 
basic  state  downstream  velocity  shear, 

'  O 

the  rate  of  working  by  wave  autocorrelation  Reynolds  stress  against 
basic  state  cross-stream  current  gradient. 


E  do  G(0)<'(u'2  + 
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the  net  advection  of  wave  kinetic  energy  by  basic  state  cross-stream 
volume  flux. 


and 


00 


the  rate  of  working  by  wave  pressure  gradient  forces.  The  angle  brackets 

here  denote  quantities  averaged  over  a  wave  period.  The  terms  on  the 

right  of  (31)  may  be  computed  by  taking  the  real  parts  of  the  products 

of  the  fluctuating  properties  involved  and  using  the  wave  forms  and 

2 

subsequent  expansions  in  k  for  the  eigenfunctions.  One  then  finds  to 
2 

lowest  order  in  k  and  F: 


''xy  “  SeF  *  V.-2  log  d.''  -  v.'’  - 

"xx  ' 

“  SgFY/2 , 

Wp  =  -SgFY(vo"^  log  do"^  -  Vo"b. 


where  y  =  exp  (2FKy)/2,  the  amplitude  growth  factor.  Thus,  in  the  net 

kinetic  energy  production  expressed  by  the  right  side  of  (31),  is 

negligible  while  the  other  three  terms  are  0(F).  R  has  the  sign  of  S 

xy  e 

for  0<Vo<l  so  that  this  Reynolds  stress  term  extracts  energy  from  the 
basic  state  flow  for  5^=1  but  supplies  it  for  Sg=-1.  The  log  singularity 
for  do->0  reflects  the  singularity  in  u(0)  for  that  limit.  The  advection 
of  wave  kinetic  energy  likewise  has  the  sign  of  Sg,  but  the  pressure 
gradient  work  term  has  the  opposite  sign.  In  summing  tnese  terms  the 
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two  log  singularity  terms  cancel  and  the  net  result  is  simply  the  first 


term  of  R. 


Consequently,  wave  kinetic  energy  growth  or  instability  occurs  with 
Sg=l ,  and  decay  with  Sg=-1 ,  consistent  with  the  result  (30)  for  wave 
amplitude  growth  k.  For  the  limiting  case  Vo-»-l  (do-^),  one  has  simply 


dK'/dt  =  S^Fy/6  =  Fky/12. 


In  the  net  production  or  decay  of  wave  kinetic  energy,  then,  the 
wave  cross-correlation  Reynolds  stress  ^-u'v'^  plays  the  critical 
role,  as  one  would  expect  from  a  model  with  only  one  active  layer 
("barotropic").  To  0(F)  and  lowest  order  in  k  one  may  write  the 
fluctuation  product  as 


-u'v'  =  -2y  sin  a  cos  a  +  FYo^*^^X^ ^^^cos^a) 


where  a  =  ky-ojt.  The  time  mean  of  the  first  term  on  the  right  vanishes 
because  its  harmonic  functions  are  in  quadrature  phase,  i.e.,  the  quadrature 
phase  between  u'  and  v'  in  the  F=0  flow  can  produce  no  Reynolds  stress. 
However,  the  0(F)  flow  alters  this  phase,  as  the  0(F)  terms  in  the  wave 
forms  (10  a,b)  for  u'  and  v'  show.  In  consequence,  the  second  term 
above  produces  the  time  mean  Reynolds  stress  of  0(F)  that,  in  turn, 
accounts  for  the  net  production  or  loss  of  wave  kinetic  energy. 


In  this  paper  I  have  developed  a  simple,  analytic  model  of  small 
amplitude,  long,  free  wave  disturbances  to  an  established  upper  ocean 
density  front.  The  major  characteristic  of  the  model  is  the  structure 
of  the  basic  state  assumed,  essentially  Stommel's  (1966)  model  of  the 
Gulf  Stream  in  deep  water  where  the  active,  upper  layer  has  uniform 
potential  vorticity.  Two  major  categories  of  frontal  wave  dynamics 
result;  the  slow,  downstream  propagation  of  long  waves  which  are  stable 
in  the  absence  of  basic  state  cross-stream  flow  (F=0),  and  modified 
forms  of  these  which  grow  or  decay  downstream  in  the  presence  of  small 
basic  state  cross-flow  (0<F<<1).  The  model  fits  within  the  general 
class  of  barotropic  models,  as  it  has  only  one  active  layer,  but  falls 
outside  the  class  of  quasi -geostrophic  models,  as  the  basic  state 
current  has  a  Rossby  number  of  0(1).  It  may  be  viewed  as  an  extension 
of  G83  for  the  long  wave  or  zeroth  mode  obtained  by  relaxing  the  condition 
of  no  shear  across  the  inner  boundary. 

The  major  properties  of  the  waves  for  F=0  are:  (1)  phase  and  group 

speeds  which  are  much  slower  than  the  basic  state  downstream  current, 

(2)  phase  speeds  which  increase  with  the  wavenumber  squared  or  with  the 

frequency  to  the  two-thirds  power,  and  (3)  anomalous  dispersion.  To 
2 

lowest  order  in  k  the  cross-stream  momentum  balance  for  the  waves  is 
geostrophic  and  potential  vorticity  changes  are  surpressed  by  a  simple 
balance  between  perturbations  in  lateral  shear  of  the  downstream  current 
and  interface  depth.  In  the  wave  troughs,  where  the  current  is  faster, 
the  interface  shoals,  and  the  reverse  obtains  for  the  crests.  The 
downstream  momentum  balance,  in  contrast,  is  neither  simple  nor  geostrophic. 


2 


I 


but  contains  both  downstream  and  cross-stream  advection  of  momentum,  as 
well  as  the  geostrophic  terms.  For  the  limiting  case  of  no  dissipative 

zone  (do-^)  the  lowest  order  results  agree  with  those  of  Paldor  (1982) 

2 

for  k  -*0.  They  are  also  consistent  with  the  numerical  model  of  Martin 
(1981).  In  application  to  Gulf  Stream  meanders  the  model  gives  phase 
speeds  or  wavelengths  that  agree  over  the  long  wave  band  with  the  results 
of  Halliwell  and  Wooers  (1983)  based  on  satellite  observations  of  propagating 
waves.  Thus,  even  for  as  deep  a  flow  as  the  Gulf  Stream,  one  need  not 
necessarily  resort  to  more  complicated  two-layer  or  baroclinic  models, 
such  as  Orlanski's  (1969),  to  obtain  similar  wave  properties  to  those 
observed. 

A  simple,  barotropic  long  wave  instability  mechanism  was  found  upon 
adding  small  cross-stream  flow  u  to  the  basic  state,  such  as  would  be 
driven  by  turbulent  mass  entrainment  across  the  frontal  interface  within 
the  narrow  dissipative  zone.  The  mechanism  differs  fundamentally  from 
that  of  classical  barotropic  instability  in  quasi -geostrophic  models 
which  require  an  extremum  in  basic  state  potential  vorticity.  Here,  in 
contrast,  the  potential  vorticity  is  uniform.  The  model  shows  that  wave 
amplitude  and  kinetic  energy  grow  for  cross-flow  toward  the  parent  pool, 
decay  for  the  reverse,  and  are  stable  for  zero  cross-flow.  The  spatial 
rate  of  downstream  amplitude  growth  is  independent  of  wavenumber  for 
long  waves  and  has  a  downstream  scale  of  0(x/F),  many  Rossby  radii  or 
frontal  zone  widths.  Finite  cross-flow  eliminates  the  perfect  quadrature 
phase  found  between  u'  and  v'  for  F=0  to  yield  small  cross-correlation 
Reynolds  stress  ^-u'v'^  which  works  against  the  basic  state  shear  of 
downstream  velocity  dv/dx  to  change  wave  kinetic  energy.  The  results 
are  consistent  with  those  of  Halliwell  and  Mooers  (1983),  especially  in 


the  lack  of  dependence  of  spatial  growth  rate  on  wavenumber. 

The  primary  object  of  the  model  is  to  illuminate  basic  frontal  wave 
mechanisms  and  to  suggest  approaches  for  more  sophisticated  models, 
especially  numerical  ones.  The  latter,  in  particular,  should  permit 
full  coupling  of  a  dissipative  inner  zone  structure,  such  as  a  time- 
dependent  version  of  Assaf's  (1977)  model,  with  a  nearly  inviscid  outer 
zone.  In  addition,  they  might  treat  other  potentially  important  compli¬ 
cations,  necessarily  omitted  here,  such  as  finite  amplitude  waves  and 
downstream  variation  in  the  basic  state,  simulating,  for  example, 
progressive  departure  of  volume  flux  from  the  Gulf  Stream  into  the 
interior,  recirculation  region. 
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